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Abstract 

In the framework of the recently developed asymptotic models for tibio-femoral 
contact incorporating frictionless elliptical contact interaction between thin elastic, 
viscoelastic, or biphasic cartilage layers, we apply an asymptotic modeling approach 
for analytical evaluating the sensitivity of crucial parameters in joint contact me- 
chanics due to small variations in the thicknesses of the contacting cartilage layers. 
The four term asymptotic expansion for the normal displacement at the contact 
surface is explicitly derived, which recovers the corresponding solution obtained 
previously for the 2D case in the compressible case. It was found that to minimize 
the influence of the cartilage thickness non-uniformity on the force-displacement 
relationship, the effective thicknesses of articular layers should be determined from 
a special optimization criterion. 

Key words: Contact problem, thin elastic layer, variable thickness, asymptotic 
model, articular contact 



1 Introduction 

Contact problems involving transmission of forces across biological joints are of consid- 
erable practical importance and a number of numerical models for articular contact are 
available [T51I22] . At the same time, the necessity of analytical models becomes an im- 
portant issue in developing improved understanding of load distribution in the normal 
and pathological joints, which affects the mechanical aspects of osteoarthritis [91123] . Also, 
analytical modelling of the distributed internal forces generated by articular contact in 
tibio-femoral joints is required in multibody dynamic simulations of physical exercise of 
a human skeleton [HJTH]. As a rule, analytical models of articular contact assume rigid 
bones and represent cartilage as a thin elastic layer of constant thickness resisting to 
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deformation like a Winkler foundation consisting of a series of discrete springs with con- 
stant length and stiffness [12] • However, a subject-specific approach to articular contact 
mechanics requires developing patient-specific models for accurate predictions. Recently, 
a sensitivity analysis of finite element models of hip cartilage mechanics with respect to 
varying degrees of simplified geometry was performed in [3]. 

Based on the asymptotic analysis of the frictionless contact problem for a thin elastic 
layer bonded to a rigid substrate in the thin-layer limit piTTO] , the following asymptotic 
model for contact interaction of two thin incompressible layers was established [5] : 

-(E^hl + E 2 1 hl)A yP (y) = 5 -ip(y), y£u, (1) 

p(y) = o, |£(y) = o, yer. (2) 

Here, p(y) is the contact pressure density, h a and E a are the thickness and elastic mod- 
ulus of the layer material, respectively, a — 1, 2, A y — d 2 /dy\ + d 2 /dy 2 is the Laplace 
differential operator, <5 is the vertical approach of the rigid substrates, <f(y) is the gap 
function defined as the distance between the layer surfaces in the vertical direction, uj is 
the contact area, V is the contour of u, d/dn is the normal derivative. 

It was shown piUfTT] that the problem ([T|, ^ describes the instantaneous response of 
thin biphasic layers to dynamic and impact loading. In [7], the elastic model Q, ^ was 
generalized for the general viscoelastic case. 

With respect to articular contact, a special interest represents the case when the subchon- 
dral bones are shaped as an elliptic paraboloid 

cp(y) = (2i? 1 )- 1 y 1 2 + (2R 2 )- l y 2 2 (3) 

with positive curvature radii R\ and R 2 . 

In the case ([3]), the exact solution to the problem ([I]), (|2| has the following form pfTU] : 

(4) 

Integration of the pressure distribution (|4]) over the elliptical contact region to with the 
semi-axes a\ and a 2 results in the following force-displacement relationship [5]: 

7TT7? 

P = —M P (s)R 1 R 2 5l (5) 

Here, Mp(s) is a dimensionless factor depending on the aspect ratio s = a 2 /ai, and the 
coefficient m is given by 

m={E^h\ + E^hl)-\ (6) 



The asymptotic model (JT[) — (JSJ) assumes that the cartilage layers have constant thicknesses, 
whereas it is well known p] that articular cartilage has a variable thickness as well as the 
surface of subchondral bone deviates from the ellipsoid shape [20]. A sensitivity of the 
model Q, (|2| with respect to small perturbations of the gap function ^ was performed 
in [5]. In particular, it was shown [3] that the influence of the gap function variation on the 



force-displacement relationship will be negligible if the effective geometrical characteristics 
Ri and R 2 are determined by a least square method. 

To our knowledge, in the literature there is only one study [21] where the 2D case of contact 
problem for a thin elastic strip of variable thickness was solved by an asymptotic method 
under the assumption that Poisson's ratio of the strip material is not very close to 0.5. 
At the same time, many asymptotic solutions were derived for an elastic layer of constant 
thickness both in the axisymmetric [2l ll3|19] and the non-axisymmetric |^|10jll6|17j cases. 

In the present paper, a three-dimensional unilateral contact problem for a thin elastic 
layer of variable thickness bonded to a rigid substrate is considered. Two cases are stud- 
ied separately: (a) Poisson's ratio of the layer material is not very close to 0.5; (b) the layer 
material is incompressible with Poisson's ratio of 0.5. After developing a refined asymp- 
totic model, we apply sensitivity analysis to determine how "sensitive" is the mathematical 
model Q, (]2]) to variations in the values of the layer thicknesses. To be more precise we 
consider the term "sensitivity" in a broad sense by allowing variable layer thicknesses, 
whereas the original model deals with scalar parameters h\ and hi. 



2 Contact problem for a thin elastic layer with variable thickness 



We consider a homogeneous, isotropic, linearly elastic layer with a plane contact surface, 
X3 = 0, and a variable thickness, H(x\,X2), firmly attached to an uneven rigid surface 

x 3 = H(x 1 ,x 2 ). (7) 



In the absence of body forces, equations governing small deformations of the layer are 

dau daoj da^j , 
^ + 7^ + 7^=0, J = 1,2,3, 

OX\ OX2 OX3 



dij = X5ij(e n + e 2 2 + £33) + 2/ie 

1 / diij du-j 
+ 



2 \dxj dxi 



(9) 
(10) 



where <Jij is the Cauchy stress tensor, A and /1 are Lame parameters for the layer material, 
Eij is the infinitesimal strain tensor, Uj is the displacement component along the Xj-axis, 
Sij is Kronecker's delta. 




*i = ;yi 



x 3 = z 



Figure 1. Elastic layer with a variable thickness. 



We assume that the elastic layer is indented by a smooth rigid punch in the form of an 
elliptic paraboloid 

x 3 = -<p(x 1 ,x 2 ), <p(x 1 ,x 2 ) = {2R l )- l x\ + {2R 2 )- 1 x 2 2 . (11) 

Under the assumption of frictionless contact, we have 

<r 3 i(y,0) = 0, a 32 (y,0) = 0, y = ( yi ,y 2 ) ER 2 . (12) 

Here the Cartesian coordinate system (yx, y 2 , z) is used such that y\ = x\, y 2 = x 2 , z = x 3 . 

Denoting by 5q the indenter's displacement, we formulate the boundary condition on the 
contact surface as follows: 

u 3 (y,0) >S - V (y), ( T 33 (y,0)<0, 
(u 3 (y,0)-6 + <p{y))a 33 {y,0) = 0, yGl 2 

On the rigid substrate surface ([7]), we will have 

Uj (y,H(y)) = 0, yGl 2 (j = 1,2,3). (14) 

Assuming that the layer is relatively thin in comparison to the characteristic dimensions 
of u>, we introduce a small dimensionless parameter e and set 

S = eS* , i?i = e- l R\, R 2 = e^R* 2 , (15) 

H(y)=eK(l + e^(y)), (16) 

where 5q, i?*, and R 2 are assumed to be comparable with h*, and, moreover, |^*(y)| < h* 
for y e R 2 . 

The problem is to calculate the contact pressure distribution 

p(xi,x 2 ) = -a 33 (xi,x 2 ,0), (xi, x 2 ) G u. (17) 

An important characteristic of the problem is the contact force required to indent the 
layer 

P=[[p(y)dy. (18) 



Finally, let us introduce the notation 

h = eh, H(y)=e 2 H*(y), (19) 

where 

#*(y) = M*(y)- (20) 



Thus, the following relation takes place: 

H(y) = h + H(y). 



(21) 



Here, h is an average thickness, H(y) is a small variation such that H(y) <C h (see 
Figure [I]) . 



3 Asymptotic Ansatz 



First, we introduce the so-called stretched coordinate 

C = e~ l z. 



(22) 



Now, substituting ^ and (10) into Eqs. Q and taking into account (22), we arrive at 
the following Lame system for the displacement vector u = (v, w): 

duo 

£ ~ 2 ^de + £ ~ 1(A + ^ Vy dc + fiVy ' v " v + (A + ^ VyVy ' v = °' (23) 



£~ 2 (2/i + ^)^+ e 0* + v)Vy ■ ^7 + l&v ■ Vj,w = 0. 



(24) 



Here, V y = (d/dy 1,8/81/2) is the nabla differential operator, and the dot denotes the 
scalar product, so that V y ■ V y = A y is the Laplace operator. 



Correspondingly, the boundary condition (12) takes the form 

_i<9v 



0. 



(25) 



In view of (13) and (15), we will have 

w(y,0) =e(5* -ip*(y)), y e u, 



where we introduced the notation (see Eqs. (11) and (15)) 

^(y) = (2 J R*)- 1 2/l 2 + (2^)- 1 l/ 2 2 . 



(26) 



(27) 



Furthermore, by stretching the vertical coordinate, formula (17) is transformed as 



p(y) =e- 1 (2fi + X)— + AV y -v 



C=o 



y Gu. 



(28) 



Finally, the boundary conditions (14) on the substrate surface (see Eqs. (|7|) and (16)) 

C = K(l + e^(y)) (29) 



take the following form: 

v(y, h* + £/i*^*(y)) = °7 w(y, K + e/i*^*(y)) = 0. 



(30) 



Observe that among Eqs. (23) -(26), (28), and (30), there are only two inhomogeneous 



ones, namely, (26) and (28). So, the form of (26) suggests the asymptotic expansion 



(31) 



Now, in view of (28) and (31), we may suggest that p(y) = 0(1) as e — > 0. However, 



taking into account the homogeneous conditions (25) and (30), we put 



v(y,0 = ev (y,() + £ V(y,() + 



(32) 



We emphasize that the asymptotic Ansatz (31), (32) is valid only inside the contact 



region u. In other words, a plane boundary layer should be constructed near the edge of 
the contact area. We refer to [jfT3] for more details. 



4 Derivation of asymptotic expansions 



Substitution of (31) and (32) into Eqs. (23) and (24) gives 



( ' dw 1 <9 2 v 2 \ 

+ e° UA^v + (A + (i) V y V y • v° + (A + ^)V„^r + (i^j + ■ ■ . = 0, 

e- 2 (2^ + A) — + e- 1 + A) — + (A + (i)V y ■ —J 
+ e° 1(2(1 + A) — + (A + (i) V y ■ — + /iA^° 1 + . . . = 0. 



(33) 



(34) 



Further, the substitution of (31) and (32) into the boundary conditions (25) and (28) at 
the contact region yields 



idv „/ n dv 1 ^ / . <9v 2 



o. 



C=o 



(2/x + A) ^ + e (\Vy • v° + (2(1 + A) 



+ £ 2 (\VyV 1 + (2(i + \)—) + 



5( 



C=o 



-p(y)- 



(35) 



(36) 



Finally, the substitution of (|31|) and (|32|) into (|30|) leads to the boundary conditions 

0. 



+ £ 2 V 2 + #, 



~ av 1 h 2 a 2 v° 



9( + 2 ac 2 



+ . . . 



(37) 



ii T ~ T dw°\ 2 / 2 ~ dw 1 H 2 d 2 w° 



c=/u 



(a 



where the notation (20) was taken into account. 



Thus, on the basis of Eqs. (33) - (38), we arrive at a recurrence system of boundary- value 
problems for the functions and w k (k = 0, 1, . . .). In the next two sections, we will 



construct the first several terms of the asymptotic series (31) and (32). 



5 Asymptotic model for a compressible elastic layer 



According to (|33|) - (|38|), the first-order problem takes the form 

dw° 



d 2 w ° 

^c 2 



(2// + A) 



o, Ce(OA), (2/i + A) 



C=o 



p(y), w 



<9 2 v° / t \ <9v° 



<9C 2 



o, v° 



C=o 



t=h* 



0. 



From (39) and (40), it immediately follows that 



v°(y,C) = o. 



(39) 
(40) 



(41) 
(42) 



In view of (42 ), the second-order problem, derived from Eqs. (33 ) - (38 ), takes the following 
form: 



d 2 w 1 
IK 2 



o, c e (o, ^ 



dw 1 



0, w 1 



C=o 



-#*(y) 



d 2 \ l dw° Sv 1 



9C 2 



C=o 



-V^w 



C=o' 



C=/i* 



(43) 



0. (44) 



It can be easily shown that the solution of the problem (43) reads as 



(45) 



On the other hand, the unique solution of the problem (44) can be represented as 

v 1 (y,C) = *(C)V^(y), 
where we introduced the notation 



*(C) 



A " -(^-C 2 ) + 7r^T(^-C). 



2/i(2/i + A) 



2/i + A 



(46) 



(47) 



tion for the normal displacement 



Thus, collecting Eqs. (31), (41), and (45), we obtain the two-term asymptotic approxima- 

(48) 



p(y) 



p(y) 



H.(y). 



By taking into account the scaling relations (19), we rewrite (48) in the form 

p(y) 



2/i + A 



(49) 



Now, substituting the expression (49) into the contact condition 



«3(y,o) = 5 - (p(y), yew, 

we derive the following equation for the contact pressure density: 

h + H(y) 



(50) 



2/i + A 



-p(y) = £o-^(y), yew. 



(51) 



In view of the condition p(y) > for y € w, we get 

2/i + A 



p(y) 



(*o - <p(y) 



h + H(y) 

where (x) + = max{x,0} is the positive-part function. 



(52) 



Now, invoking the notation (|21|) for the variable thickness of the elastic layer, we rewrite 

2/i + A 



(52) as follows: 



p(y) 



H(y) 



(5o - <p(y))_ 



(53) 



Formula (53) shows that a thin compressible elastic layer deforms like a Winkler founda- 



tion with the variable foundation modulus 



My) 



2/i + A 



(54) 



Finally, let us recall that the Lame parameters A and \i are related to Young's modulus, 
E, and Poisson's ratio, v, by formulas 



A 



Ev 



E 



;i + 2i/)' ^ 2(1 + 1/)" 



(55) 



Now, in view of (55), it is readily seen from (54) that k(y) — > oo as v — > 0.5. This 
implies that the case of an incompressible elastic layer with v = 0.5 requires a special 
consideration. 



6 Asymptotic model for an incompressible elastic layer 



Let us continue the process of constructing terms of the asymptotic expansions (31) and 



(32). In view of (45), Eqs. (33) -(38) yield the third-order problem 



2 w 2 0~v^ 
(2/j + A)^- = -(A + fi)V y ■ ^— - fiAyiu , C e (0, h* 



(2fi + A) 



dw 2 



<9 2 v 2 



o, Ce(o,U 



«9v 2 



-AV„ ■ v 1 



(=0 



C=o 



C=o' 



w 



-VyW 1 



C=0' 



C=h t 



dv 1 



(56) 
(57) 



Substituting (41) and (46) into Eqs. (56), we derive the boundary-value problem 

d 2 w 2 



d( 2 
dw 2 



AA y p(y) 
fi(2fi + A) 2 

X(ji - X)hl 



C=o 



2//(2/x + A) 



[(2fi + X)C-tiK], £e(o,K) 

-A y p(y), w 2 



C=h* 



0. 



(5f 



It can be checked that the solution to ( 58 ) can be condensed to the form 



w 



- A(2/i + A)(C 3 - hi) - 3\(ji - X)hl(C - K)}. 



(59) 



Further, in view of (45) and (46), the problem (57) takes the form 
«9 2 v 2 



d( 2 



o, (6(o,4 



«9v 2 



Xh» 



<=h. fi(2fi + A) 



C=o 2/i + A 

Here the arguments of functions p(y) and if* (y) are omitted for clarity. 
It can be easily verified that the solution to ((601) has the form 



H.V v p. (60) 



v 2 (y,C) = £rAv v (pH.) - , Afe * ^ H*V yP . 



2/i + A 



/i(2// + A) 



(61) 



We emphasize that in contrast to the first two term approximation (48), the third term 



(59) does not vanish at the contact surface in the limit as v — > 0.5. Indeed, formula (59) 



yields 



2/ m ^A(A-/i) 
w(y '° ) = ~ 3M2/i + A) 2 A ^ (y) ' 



(62) 



where, in view of (55), as v — > 0.5, we will have 



X(X-H) v{l + v){Av-l) 3 



/i(2/i + A) ; 



£7" 



(63) 



Finally, in order to construct a correction for the leading asymptotic term (59 ), we consider 
the following problem: 

d 2 w 3 d~v 2 
(2/^ + A) — — = -(A + /x)Vj,- -/AX> (0,/i*), 



<9C 2 



(2/i + A) 



<9-ur 



<9C 



C=o 



-AV„ ■ v 2 



C=o' 



w 



(=h t 



-H*(y) 



dw 2 



C=h, 



(64) 



Substituting the expressions (45), (59), and (61) into Eqs. (64), we get 

d 2 w 3 



A 



9(2 - ( )a A(pff.), Ce(o,K) 



dw 3 



C=o 



A/i* 
At(2/x + A) 



(AV a ■ (£>*V y p) -nA y (pH*)) , 



10 



Xhl 

t=h* 2(2/2 + A) 2 



(65) 

#*Aj,p. (66) 



Integrating Eq. (65), we obtain 



w 3 (y,0) 



A 



A(p^)C 2 + ^i(y)C + Co(y), 



(67) 



2(2/i + A) 2 yv 

where the integration functions Ci(y) and Co(y) are restricted by the boundary conditions 



(66). So, it can be checked that 



Ci(y) 



A/i* 



M 2 /U + A) 



XV y • (H*V yP ) - fxA y ( P H, 



C (y) 



A/) 2 ~ \ 2 h 2 



2(2/x + A) 2 ' 



/i(2/i + A) 



From (67), it immediately follows that 

W 3 (y,0) = C (y). 



(68) 
(69) 

(70) 



Note that the elastic constants entering Eq. (69) are evaluated in terms of the engineering 
elastic constants as 



A v(l + v)(l-2v) A 2 



(2^ + xy 



E(l - u) 



2v 2 a + v) 



/i(2/i + A) 2 



,2 ' 



(71) 



Hence, in view of (71), the first term in (69) disappears as v — > 0.5, and we get 

3h 2 



C (y) 



v=0.5 



£7 



(72) 



Thus, collecting Eqs. (31), (62), (70), and (72), we obtain the following two-term asymp 



totic approximation for the normal displacement at the contact surface in the case of 
incompressible elastic layer: 



h 3 3h 2 
MY, 0) * s 3 ^A y p(y) - e^V y ■ (H*(y)V y p(y)). 



Recollecting the scaling relations (19), we rewrite (73) in the form 

3h 2 



« 3 (y,0) * -^A y p(y) - ^V y ■ (H(y)V y p(y)). 



(73) 



(74) 



Now, substituting the expression (74) into the contact condition (50), we arrive at a 



partial differential equation in the domain u with respect to the function p(y). According 
to the asymptotic analysis [13] . at the contour T of u, we impose the following boundary 
conditions: 

p(y) = o, ^(y) = o» yer. (75) 

Here, d/dn is the normal derivative. We stress that the location of the contour T must 
be determined as part of the solution. 



7 Comparison of the obtained result with the solution for the 2D case 



Collecting Eqs. (31), (41), (45), (59), and (67), we obtain 



w{y,0)~e 



h* 



2/i + A 



P(y)^ 2 ^4p(y)- 3 ^ A(A " M) 



+ e q 



2/i + A" w/ " 3/i(2/i + A) 2 y - 
H4y)A y p{y) + A y (p(y)H4y) 



A y p(y) 



[2 (2/i + A) 2 

V/r rV, • (£*(y)V yP (y) 



/t(2/t + A) 2 



(76) 



Substituting the asymptotic expansion (76) into the contact condition (26) and using the 



notation (20), we derive the following equation for the contact pressure density: 



P(y) + #*(yMy) - e 2 ^ A , ^ A yP (y) 



3/x(2/x + A) 



+ e 



[I 



2/i(2/i + A) 
-2AV„ - (v.(y)V v p(y) y /; 



^*(y)A 2/ p(y) + A y (p(y)^(y 



2 " + A /'(y). 



Here we also introduced a shorthand notation for the right-hand side of (26), i. e. 

r(y) = *;-¥>* (y)- 



(77) 



(7* 



It should be noted that Eq. (77) is applied for the case of compressible materials when 



its right-hand side makes sense. 



By applying a perturbation method, a solution to Eq. (77) is represented in the form 



p(y) - 2/ fr + A (°"°( y ) + £(J i(y) + £2(J 2(y) + eV 3 (y; 



(79) 



After the substitution of (79) into (77), we straightforwardly obtain 



JO) 



^ = -vsr 



^A(A - /i) 
3/i(2/i + A) 



(^A,r + a,(^d) 



h?A 



2fi(2fi + A) 



where for the sake of brevity, the argument y is omitted. 
Further, making use of the differential identities 



A(/V>)=^A/ + /A^ + 2V/-V^, 



we simplify formula (81) as follows: 



as = + 3M2/X + A) 



(v,r.v^ + ^ r) -^±|rA^ 



52) 



Now it can be easily checked that the four-term asymptotic expansion (79) with the 



coefficients given by (80) and (82) in the 2D case recovers the corresponding solution 



obtained in [21], where also the next asymptotic term in (79) was explicitly written out. 



Finally, observe that formula (69) can be transformed into the following one: 



C (y) 



[V y H* ■ V y p + H*A y p\ + — — ——pAyH*. 



ti(2fi + A) 2 



2(2// + A) 2J 



(83) 



Note also that the expression in the brackets in (> 



83) is equal to V y ■ (H*V y p). 



8 Application to sensitivity analysis of articular contact mechanics 



According to (74), (75) the refined asymptotic model for contact interaction of thin in- 



compressible layers bonded to rigid substrates looks as follows: 

- m^Aypfr) - ifV y ■ (H a (y)V y p(y)) = 5 - <p(y), yea), 

a=l ha 

p(y) = o, ^(y) = o. yef. 

Here, f is the contour of the contact region a), and the notation (jfjl) was used. 
Let us put 

p(y) =p(y) +p(y), 

where p(y) is the solution to the original asymptotic model (fTl), pj. 



(84) 
(85) 



36) 



Then, under the assumption that the thickness variation functions H\(y) and -f^(y) 
introduce a small variation into the elliptical contact region u corresponding to the density 



p(y), we derive from (84) - (86) the following limit problem for the variation of the contact 
pressure density: 

2 3/> 2 

- m-'Aypiy) = £ _H Vy ■ (F a (y)V^(y)), y € u, (87) 

p(y) = 0, yer. (88) 
Here, T is the contour corresponding to the contact pressure Q. 

Moreover, the thickness variations Hi(y) and -f^(y) will not greatly influence the resulting 
force-displacement relationship, if 



JJp(y)dy = o. 



39) 



Let us derive the conditions for Hi(y) and H 2 (y) under which the equality (89) holds 
true. With this aim we consider an auxiliary problem 



A,0(y) = l, yeu, 6(y) = 0, yeT 
with the solution 0(y) = — a\a\{2{a\ + al)y 1 0{y), where 



0(y) = l 



a\ a 2 . 



(90) 



(91) 



In view of (90), we rewrite Eq. (89) as 



JJp(y)A y e(y)dy = 0. 



(92) 



Now, applying the second Green's formula and taking into account Eqs. (87), (88), and 
(90), we reduce Eq. (92) to the following one: 



W E ■ (H a (y)V y p(y)) dy = 0. 



(93) 



After rewriting Eq. (93) in the form 

2 



d ( dp 



0(y) E n , ., 



(y)E^«(y)W = o 

a=l / 



and integrating by parts with (91 ) taken into account, we get 



1 2^ dp_ 
3=1 a l <% 



2 /i 2 



+ / 6(y) cos(n,^)|f (y) £ f^ a (y) ds y = 0. 



(94) 



It is clear that the line integral in (94) vanishes due to the boundary condition (91 ). Hence 



taking into account the exact expression Q for p(y), we finally transform Eq. (94) into 
the following one: 

2 hi rr ~ 

(95) 



Eir II H a (y)p(y)dy = 0. 

a=l 



Here we introduced the notation 



p(y) 



'syt + yV\ ( l _yl_yJ 

a\ sa\)\ a\ a\ j 



(96) 



Based on the derived Eq. (95), we suggest the following optimization criterion for deter- 



mining the average thicknesses h\ and /i2 : 

min / / (H a (y) - h a ) 2 p*(y) dy. 



(97) 



Here, is an characteristic elliptic domain with semi-axes a* and b*. In particular, in 
the capacity of one can take the average contact area for a class of admissible contact 
loadings, and p*(y) is given by 



p*(y) 



' s*yf 

of ^ s*al 2 



1 - 



vl vl 



*2 i 



(9£ 



where s* = a\ja\ is the aspect ratio of w*. 



It is clear that the necessary optimality condition for (97) has the form 



fj(H a (y)-K)p*(y) dy = 0, 



(99) 



from where it follows that 

ha = [II p *^ dy ) II H *(y)p*(y) d y- 



(100) 



It is left to show that Eq. (]95|) follows from (|99|) if Co>* coincides with u>. Indeed, in view of 

(101) 



(21), Eq. (99) is equivalent to the following one: 

jj H a (y)p*(y)dy = 0, 



a = 1,2. 



Now, adding the two equations above multiplied by E a h a) a = 1,2, respectively, we 



arrive at Eq. (95). 



It is interesting to observe that Eq. (100) indicates that in order to obtain the optimal 



average thickness h a , the corresponding variable thickness H a (y) has been averaged with 



the weight function p*(y) given by (98). 



Finally, note that in the case of compressible layers, the optimal value of the average 
thickness h a coincides with the simple average of H a (y). 



9 Discussion 



Let us compare the optimal condition (101) with the condition derived in j5] for the 
problem of gap variation, for which the limit problem for the variation of the contact 
pressure density takes the form 



m 1 A y p(y) = <p(y), yew, p(y) = 0, y G T, 
where y5(y) is the ellipsoidal gap variation. 



(102) 



In this case, the optimization condition (89), which is conveniently rewritten as 



Hp(y)A y 6(y)dy = 0, 

U! 

after applying the second Green's formula can be replaced with the following one: 

II v(y)0(y)dy = 0. 



(103) 



Finally, with the reference to the characteristic contact region a;*, Eq. (103) takes the 
form 

// <p(y)e*(y)dy = 0, (104) 



where we introduced the notation 



0*(y) = i 



y±_ 

~ #2 



vl 



(105) 



Thus, comparing Eqs. (101) and (104), we see that they differ only by their weight func- 



tions p*(y) and #*(y), respectively, given by (98) and (105). It is interesting to observe 



that the maximum of the function p* (y ) does not coincide with the center of the domain 
a;*, while the function #*(y) makes an emphasis namely on the central part of u*. 

It should be noted that the least squares optimization criterion (which is analogous to 



the integral (97) without the weight function) derived in [5] from an ill-posed problem 
underestimates the contribution of the central part of the contact domain u. Thus, to 
minimize the effect of the peripheral part of u, it is required to decrease the characteristic 
contact region u*. 



10 Conclusion 



The present study results in asymptotic solutions to the three-dimensional unilateral 
contact problem for a bonded thin elastic compressible or incompressible layer of variable 
thickness. As the main result of the present paper, the four term asymptotic expansion 



(76) is obtained for the normal displacement at the contact surface. In the compressible 



case, the derived asymptotic expansion (79) for the contact pressure state recovers the 



corresponding solution obtained in [2TJ for the 2D case. 

The objective of this study was to apply an asymptotic modeling approach for evaluating 
the sensitivity of the asymptotic model of tibio-femoral contact due to small variations in 
the thicknesses of the contacting articular cartilage layers. It was found that to minimize 
the influence of the cartilage thickness non-uniformity on the force-displacement relation- 
ship, the effective geometrical characteristics hi and of articular layers, which enter 
the asymptotic model Q, (j2|, should be determined from the introduced optimization 
criterion (97). 
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